Symmetry properties of the nodal superconductor PrOs 4 Sbi 2 
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We present a theoretical study of the superconducting gap function in PrOs4Sbi2 using a 
symmetry-based approach. A three-component order parameter in the triplet channel best de- 
scribes superconductivity. The gap function is non-degenerate and the lower branch has four cusp 
nodes at unusual points of the Fermi surface, which lead to power law behaviours in the density of 
states, specific heat and nuclear spin relaxation rate. 
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I. INTRODUCTION 

By most accounts, PrOs4Sbi2 is an unconventional 
super conductor t i i 2 i S A 5ifiiLSi& The superconducting phase 
breaks time-reversal symmetry^ and the paired electrons 
are in a spin triplet configuration^ The existence of 
point nodes in the superconducting gap function is in- 
dicated by power law behaviour in the temperature de- 
pendencies of specific heatj2i£ penetration depth^ ther- 
mal conductivity 3 -, and Sb-NQR^ however other ex- 
periments find the gap function to be nodelessj 11 ' 12 ! 13 
Two distinct features in the specific heafr ^ 14 ' 15 and other 
measurement o 3 ' 5 ' 16 i 17 ' 18 i 19 were initially interpreted as 
two phase transitions involving a change in symmetry 
of the superconducting order parameter, but recently 
these results have been ascribed to sample inhomogeneity 
or two-band superconductivity ] 20 ! 2 1 1 22 1 23 On the theoret- 
ical side, several phenomenological unconventional order 
parameters have been proposed 24 ' 25 i 26 ' 27 and unconven- 
tional pairing mechanisms have been studied i 28 i 29 In light 
of all these intriguing and somewhat contradictory find- 
ings, it is not surprising that the only consensus on the 
symmetry of the superconducting order parameter is that 
it is probably unconventional. 

In this paper, we will consider the results of a strict 
analysis of symmetry and symmetry-breaking described 
by Landau theory. 30 i 31 i 32 According to this approach, the 
order parameter which describes the normal to super- 
conducting phase transition must belong to one of the 
irreducible representations of the crystallographic point 
group. Each irreducible representation yields a limited 
number of superconducting phases. The most convenient 
and accurate way to label the various phases is by their 
symmetry groups. All of the superconducting symme- 
try groups are subgroups of the normal phase symmetry 
G x U x JC, where G is the point group of the crystal, U 
is U(l) gauge (phase) symmetry and JC is time-reversal. 
Some of the subgroups include elements which are non- 
trivial combinations of phases, time reversal and point 
group elements. As described by Sigrist and Ueda 3 -^ and 
Volovik and Gor'kov) 3 ^ strong spin-orbit coupling is as- 
sumed in this classification scheme. 

The point group symmetry of PrOs4Sbi2 is Th (tetra- 
hedral), which has a one-dimensional representation 



A gM , a two-dimensional representation E g>u and a three- 
dimensional representation T g ^ u , in each of the singlet 
(subscripted by g) and triplet (subscripted by u) chan- 
nels. The A g order parameter describes a "conventional" 
or "s-wave" superconductor. It is associated with a sin- 
gle, fully gapped superconducting phase. The A u order 
parameter describes triplet superconductivity, also with 
a single, fully gapped superconducting phase. The A g<u 
phases have symmetry T x JC, where T is the tetrahe- 
dral point group. The E g>u and T g tl order parameters 
are each associated with more than one superconducting 
phases, corresponding to different symmetries. The E g ^ u 
order parameters describe three different superconduct- 
ing phases, of which two are accessible from the normal 
state via a second order phase transition, while the T QtU 
order parameters describe nine different superconduct- 
ing phases, of which four are accessible from the normal 
state. The symmetry properties of all of these states and 
their corresponding gap nodes are given in Table I of 
Ref. |32 
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The order parameter which best describes experiments 
is T u , the three component order parameter in the triplet 
channel. Broken time reversal symmetry rules out the 
A g ^ u order parameters. The E g>u phase that is accessi- 
ble from the normal phase and that breaks time rever- 
sal symmetry is T(D-2), which has point nodes in the 
(111) directions which are not indicated in any experi- 
ment. The T g phases which are accessible from the nor- 
mal state have either time reversal symmetry, line nodes, 
or nodes in the (111) directions, leaving T u as the only 
possibility. There are two T u phases accessible from the 
normal phase that break time reversal symmetry: C3 (E) 
and D2(E); the former has nodes in the (111) direc- 
tions, leaving the phase D2{E) as the most likely can- 
didate. The elements of the symmetry group D2(E) are 
{E, Cf/C, t/i(7r)C|/C, f/i(7r)C|}, where E is the identity, 
Ui(n) are phases, C\ are rotations of 7r about the i-axis, 
and JC is time reversal. The triplet D2(E) phase has four 
point nodes in the [±a, ±/3, 0] directions. The proof that 
D2(E) has nodes in the triplet channel is given in the 
Appendix. 

The issue of whether there are two different su- 
perconducting phases (as suggested by specific heat 
and thermal conductivity experiment a 3 ' 4 ' 14 i 15 ' 16 i 17 ) or 
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only one (according to the two-band superconductivity 
scenario 20 i 21 i 22 i 23 ) is to some extent by-passed by a fluke 
of Landau theory: the D 2 (E) phase is accessible via sec- 
ond order phase transitions both directly from the normal 
phase and via an intermediate phase D 2 {C 2 ) x IC. Thus 
it is a viable candidate for either situation. Therefore, 
we identify D 2 {C 2 ) x JC as the 'A-phase' and D2(E) as 
the 'B-phase', and we will consider both the case when 
the A-phase is present and the case when the A-phase is 
absent on the phase diagram. Note that the elements of 
the group D 2 {C 2 ) x K are {E, Cf , U(n)C%, U(tt)CI} x K 
and that D 2 {E) is a subgroup of D 2 (C 2 ) x IC. 

Recently, microscopic weak coupling theory has been 
applied to tetrahedral superconductors^^, and it was 
shown that the phase D 2 {C 2 ) x K, is stable, while D 2 (E) 
is noti^£ This is apparently in disagreement with the ob- 
servation of broken time reversal symmetry, which means 
either that PrOs4Sbi2 is a str ong co upling superconduc- 
tor, as claimed in Refs. ll3lll5lll9H22l or that the B-phase 
is better described as a D 2 {C 2 ) x JC phase. We shall not 
pursue this possibility here, apart from noting that there 
are still issues whose resolution may change the conclu- 
sions of this work. 



II. THE SUPERCONDUCTING GAP 
FUNCTION 

The superconducting gap function is a 2 x 2 matrix in 
pseudospin space, 



A(fc) = ia y ifj(k) = 
in the singlet channel, and 
A(fc) = i[a-d(k)]a y = 



ip(k) 
-if>(k) 



(1) 



-d x (k) + idy(k) d z (k) 

d z (k) d x {k)+id y (k) 

in the triplet channel, where ip(k) and d(k) are even and 
odd functions of k, respectively. For singlet pairing, the 
gap function is given by 

A(fc) = |V(fc)|, (3) 

while for triplet pairing the gap function may be non- 
degenerate, 

A±(fc)= [|d(fc)| 2 ±|g(fc)|] 1/2 , (4) 

where q(k) = id(k) x d*(k). When d(k) is real q(k) 
vanishes and the gaps are degenerate and unitary. Oth- 
erwise, the gap is non-degenerate and the lowest energy 
branch has a cusp where the two branches meet. 

The gap function may be expanded in terms of the 
basis functions for a single representation of the point 
group, 



■00) = ^2mM k ) 

i 

d{k) = ^riidi(k) 



(5) 
(6) 



where ipi(k) and di{k) are basis functions for even (spin- 
singlet) and odd (spin-triplet) representations of the 
point group, respectively, and r\i are components of the 
order parameter. For the remainder of this article we 
will limit our discussion to the three component order 
parameter in the triplet channel T u . An appropriate set 
of basis functions for this representation are 3 - 2 - 



d 1 ~ akyZ + bk z y, 
d 2 ~ ak z x + bk x z, 
ds ~ ak x y + bk y x. 



(7) 



where a and b are arbitrary real numbers. More general 
forms, which include higher orders in fc, are considered 
in the Appendix. 

The phases associated with each representation are 
minima of the Landau potential, which is expanded in 
terms of the order parameter. The transformation prop- 
erties of the basis functions (0 get transferred to the or- 
der parameter, and the Landau potential is constructed 
to be invariant under all operations of the space group, 
gauge transformations and time reversal. The Landau 
potential also determines which phases are connected by 
second order phase transitions. A complete analysis of 
the Landau potentials for the tetrahedral point group T 
is given in Ref. l32l The three component order param- 
eter (771,772,773), defined by © and ([7]), has four phases 
which are accessible from the normal state by a second or- 
der phase transition, (0, 0, 1), (1, 1, 1), (1, e 2 ™/ 3 , e" 2 "/ 3 ) 
and (0,i\r] 2 \, \r]i\), with symmetries D 2 (C 2 ) x IC, C3 x /C, 
Cs(E) and D 2 (E) respectively. Thus the components of 
the order parameter in the A-phase are (0, 0, 1) (or, more 
precisely, (0, 0, (771!)) and in the B-phase are (0, i\r] 2 \, |?7i|). 
These statements are summarised in Table I. Differ- 
ent domains of each phase are obtained by permuting 
the components; the analysis below uses this particular 
choice of domain. A discussion of domains appears in 
Section IVl 



phase 


normal - 


A 


-> B 


OP components 


(0,0,0) - 


- (0,0,177x1) 


-» (0,i[7?2|,M) 


symmetry group 


T h x U x IC - 


-> D 2 (C 2 ) x IC 


-> D 2 (E) 



TABLE I: Order parameter (OP) components and symmetry 
group elements for the proposed normal— A— >B second or- 
der phase transition sequence. Note that the A-phase can be 
skipped, since the B-phase is also accessible from the normal 
phase by a second order phase transition. 



The gap function (f4]) in the A-phase, 

A ± (k) = \ Vl \[aX + b 2 kl] 1/2 , 



(8) 



is unitary (degenerate) and has cusp point nodes in the 
[001] directions, as shown in Fig. la. In the B-phase, the 
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gap function is 



A±(fc) 



(\m\ 2 b 2 + \ m \ 2 a 2 )k 2 x + \ m \ 2 a 2 k 2 y + \ m \ 2 b 2 kl 



±2 



m\\V2\\k x \Ja 2 b 2 kl + a 4 fc2 + b*k 2 



1/2 



(9) 



In this case the gap function is non-unitary and degener- 
ate only where d(k) x d*(k) — 0, that is, along the line 
k x = 0. The gap has four nodes which are solutions to 
A_(fe) = 0. When |ryi| 2 6 2 > \ m \ 



at k v = and \/|?7i| 2 & 2 — 



'72 



and when 



the nodes are found 
2 a 2 k x = ±\i]2\bk Zl shown 

2 b 2 



< \V2\ or they are 



in Figs. lb)-ld) 
found at k z = and ^/^pa 2 — \rji\ 2 b 2 k x = ±\rj2\bk y . A 
three dimensional rendering of the lower branch of the 
gap function is shown in Fig. 2. 



a, 
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a) 




c) 






FIG. 2: (Colour) The gap function drawn over a spherical 
Fermi surface for the a) A-phase and b) B-phase. In a) the 
gap function is unitary and degenerate. In b) the gap function 
it is non-unitary and non-degenerate. Only the lower branch 
of the gap function A„(fc) is shown. 



e) 




FIG. 1: (Colour) The gap function A±(fc) drawn over a spher- 
ical Fermi surface (bold) in the k x -k z plane. In a) the gap 
function (red) is unitary and degenerate. In b)-e) the gap 
function it is non-unitary and non-degenerate. The lower 
branch A_(fc) (red) and the upper branch A+(fc) (blue) are 
both shown, a) A-phase, 772 = 0. b) B-phase, \rj2\a = 0.1|r/i|6. 
c) B-phase, \rj2\a = 0.5 171 1 £>. d) B-phase, |?72|a = 0.9|7?i|b. e) 
B-phase, \r12\a = | T71 1 fc>. 

As discussed in the Introduction, the B-phase may 
evolve either from the A-phase, with 1 772 1 <C |?7i|, or di- 
rectly from the normal phase, in which case 1 772 1 ~ 
We now discuss these two scenarios in detail. 

The order parameter of the A— >B transition is r]2, 
which increases continuously from zero at the phase tran- 
sition. The two degenerate cusp nodes in the [001] di- 
rections in the A-phase (Fig. la) split into four non- 
degenerate cusp nodes in the B-phase at the phase tran- 
sition (Fig. lb). 



The order parameter of the normal^B transition is 
\Vi\ = I V2 1- In this case, the B-phase resembles the D±(E) 
phase of octahedral systems corresponding to the three- 
dimensional representations with components (0, i, 1). In 
the Landau potential, the difference between octahedral 
and tetrahedral appears only in sixth order and higher 
terms in the order parameter^ 2 - Near the normal-to- 
superconducting phase transition, when all components 
of the order parameter are small, the growth of the or- 
der parameter is governed by fourth order terms in the 
Landau potential, which are identical for octahedral and 
tetrahedral systems, so 1 772 1 = at the phase transi- 
tion in both cases. The difference between the gap func- 
tions of octahedral and tetrahedral systems with 3D or- 
der parameter components (0,i, 1) is due to a difference 
in the basis functions ([7]): \a\ = \b\ in octahedral sys- 
tems. Thus the octahedral phase (0, i, 1) has two non- 
degenerate smooth nodes in the [100] directions shown 
in Fig. le), while the tetrahedral system has four cusp 
nodes (Figs. lb-Id). 

Thus the main difference between the two possible sce- 
narios is the positioning of the nodes at the onset of the 
B-phase. In the normals A^B scenario, the nodes will 
always be found in pairs near the [001] directions (Fig. 
lb), while in the normal— >B scenario, the positions of the 
four nodes are arbitrary (Fig. lb- Id) and depend on the 
parameters a and b. 
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III. DENSITY OF STATES 

The low temperature form of the density of states 
(DOS) in superconductors is governed by the presence of 
nodes i 31 i 37 ' 38 In general, cusp-like point nodes give rise 
to a quadratic dependence on energy. 

The DOS is given by2i 



(10) 



where E±(k) = \Je 2 {k) + A|(fc) and e(k) = - E F 
is the free particle energy. 



A. A-phase 

The gap function of the A-phase J5J) is unitary and non- 
degenerate (Fig. la). Since the main contributions to the 
integral come from the vicinity of the nodes, the integral 
over k can be split into two separate regions centred over 
each node, which are cut off such that the total integrated 
region in fc-space equals the Brillouin zone^ The nodes 
are degenerate and the contributions from each node are 
equal, 



Eq. [13] is the density of states of the phase (0,0, 1) 
{D^C^) x fC) at low temperatures. However, according 
to the considerations outlined in Section I, this phase is 
identified as the A-phase, which is only found in a narrow 
region of phase space just below H C 2- Therefore, Eq. [T3l 
is not expected to be observed in PrOs4Sbi2. 

B. B-Phase 

In triplet, non-unitary phases, in general, the gap func- 
tion is non-degenerate, except along some lines on the 
Fermi surface. All nodes are found in the lower energy 
branch of the gap function A_, and the higher energy 
branch A + is usually neglected. However, if the nodes 
are found near the line where the gaps are degenerate 
then both gaps should be taken into account. 

To find the density of states in the B-phase, we 
should consider the two different scenarios, normal— > A- 
phase— >B-phase or normal— >B-phase, separately. In the 
former scenario, 1 772 1 "C |?7i|, and pairs of nodes are found 
on opposite sides of the Fermi surface. The partners in 
each pair are very close to each other and close to the 
gap degeneracy line, as shown in Fig. lb). In this case, 
the higher energy gap should not be neglected. In the 
normal— >B-phase scenario, the positions of the nodes de- 
pend on the parameters a and b which are arbitrary. 



N(uj) 



4v 2 



dk±5(uj- E(k\\,k±)) 



dkuku 



(11) 



where fcy and kj_ are the momenta parallel and perpendic- 
ular to the Fermi surface at the node, v 2 k 2 = \r]i | 2 (a 2 fc 2 + 

b 2 k 2 ), k± = k z — k F and E(k\\,k±) « \Jk\v 2 F + k^v 2 . 

Changing variables again and using pi = vpk±^ — p cos#, 
P2 = v fen = p sm8, we find 



normal— + A-phase-^ B-phase 



At the onset of the B-phase 1 772 j "C |?7i|, and we will 



assume that \t}\\ b 2 > 



\m\ 2 a 2 . 



Then the pairs of nodes 



are found in the vicinity of [001] in the plane k y 



0, 



as shown in Fig. lb). The integration over fc-space is 
divided into four regions, which overlap for nodes within 
a pair. 

The gap function in the vicinity of the nodes for the 
case when |r/i| 2 & 2 > 1 772 1 2 a 2 can be approximated by 



N(ui) 



ab\r n \ 2 (2Tr) 2 v F Jo 

1 / 

ab\in\ 2 {2^) 2 v F J 

2co 2 



ab\rji \ 2 tt 2 vf 



dp 2 P2 



sm( 



dp\8{u - E(pi,p 2 )) 



A(fe) w V\m\ 2 b 2 -\V2\ 2 a 2 ^Jk' 2 + k' 2 (14) 



„\ ft ^where k' y = f k y and 



dpp 2 8{uj-p) (12) 



(13) 



k'u = 



V\m\ 2 b 2 - \m\ 2 a 2 
\m\b 



\m\a- 
\m\b 



(15) 



where the cutoff pq is finally introduced in the last equa- 
tion. This result is equivalent to the usual result for a 
degenerate cusp node, N(oj) = w 2 /-k 2 vfv 2 ^- apart from 
a factor of two because there are two degenerate nodes in 
our calculation. In our case, the gap velocity, defined by 
v g = VfeA(fc) is not the same in all directions since the 
node is not rotationally symmetric, and so the geometric 



average v g 



\ m \(ab) 



1/2 



appears. 



With this approximation, the '— ' branch of the gap func- 
tion continues smoothly to the '+' branch of the gap 
function at the line where the gap function is degenerate. 
Then two difficulties are overcome at once: both branches 
of the gap function are taken into account, and the con- 
tributions from each integration region are distinct, even 
though the regions overlap. Each region yields the same 
contribution to the density of states, 
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N(u) 



(2 7 r) 3 a|77 1 | 2 6 2 - |?7 2 | 2 a 2 



2 71 



dk±S(u-E(k n ,k±)) 



r 



(16) 



where v 2 k 2 = (\ m \ 2 b 2 - 1 772 1 2 a 2 ) (A;^ + ^ 
fci = Vl " l|2 ,t 2 7J' ?2|2a2 fc z T l j^k x and S(fc||,fc x ) « 



k 2 ^v 2 as before. Then performing the same 

change of variables as in the A-phase calculation, we 
find 



N{w) 



2co 2 



ir 2 VF{\r]i\ 2 b 2 — \r]2\ 2 a 2 ) 



(17) 



Note that in the limit \r/2\ — > we recover the A-phase 
result, as expected. 



2. normal — > B-phase 

In this situation, near the phase transition we have 
l^i I ~ 1*72!) however the positions of the nodes depend 
on the parameters a and b, which are completely unde- 
termined. Then there are three possibilities to consider. 
The first is shown in Fig. lb), where the nodes appear 
in pairs such that the pairs are close to the gap degener- 
acy line (if \a\ <C |6| or \b\ <C \a\); in this case the above 
calculation is valid and the result (fl~7|) is obtained for 

\m\ ~ \m\, 



N(uj) 



2a, 2 



a Tr 2 VF\rji\ 2 {b 2 — a 2 ) 



(18) 



Second, when all four nodes are spaced far apart as shown 
in Fig. lc), then the above calculations are again valid 
and the result (fT5| is obtained. 

Finally, the nodes may appear in pairs which are far 
away from the gap degeneracy line, as shown in Fig. Id). 
In this case the above treatment is invalid. Here we have 



a crossover between N(lu) ~ to 2 and N(u>) 



which is 



the behaviour of the limiting case shown in Fig. le), i.e., 
the octahedral phase (0, i, 1), with smooth (quadratic) 
nodes. Such behaviour is not observed in experiments, 
which could mean that either the components of the or- 
der parameter are unequal (normal— > A-phase— >B-phase 
scenario) or a ^ b. 



IV. SPECIFIC HEAT AND NUCLEAR SPIN 
RELAXATION RATE 

The specific heat at low temperatures is given by 3 * 



2 f c 

C[T) = tL 



duJLu 2 N(uj) 



M 

dui 



(19) 



Eqs.CiJ and Q7| yield 

C(T): 

for the A-phase, and 
C(T) 



14 7T 2 



l5vFab\i]i I 



14 tt 2 



rp6 



a 15vjr(\rii \ 2 b 2 — \rj2\ 2 CL 2 ) 



(20) 



(21) 



for the B-phase. 

The longitudinal nuclear spin-lattice relaxation rate is 
given by^i 



1)T C 



2— / duN(u)N(uj-uj Q ) 
T c Jo 



8uj 



(22) 



In the limit of small nuclear resonance frequency lu 0i one 
finds 



(1/^ 



IT 



28 



(1/Ti) Tb 15^v 2 F a 2 b 2 \ m \* T c 

in the A-phase, while in the B-phase it is 

(l/Tx) T _ b 2 28 T 5 

(l/ri) Tc ~ ^ l5^v 2 F (\ m \ 2 b 2 - \ V2 \ 2 a 2 ) 2 T^- 



(23) 



(24) 



These expressions give the low temperature behaviour 
of the specific heat and nuclear relaxation rate in terms 
of the tetrahedral parameters a and b and the order pa- 
rameter components r)i and 772- 



V. DOMAINS 

Directional dependent measurements are the ideal way 
to observe the anisotropy of the gap function. However, 
such measurements may be confounded by the presence 
of domains, different regions in space where the compo- 
nents of the order parameter are interchanged. In this 
section we offer a brief discussion of domains for the A- 
phase and the B-phase. 

The A-phase has three different domains (1,0,0), 
(0,1,0) and (0,0,1), which, in the absence of unusual 
crystal shape or external fields, are all expected to be 
present, and will lead to the observation of the full tetra- 
hedral symmetry Six (degenerate) nodes will be ob- 
served in the directions (00 ± 1). Now let us suppose 
that there is some kind of external effect along the z-axis 
which effectively lowers the symmetry from to D 2 h- In 
an octahedral system, either the single domain (0,0,1), 
with nodes in the [00 ± 1] directions will be favoured, 
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or the other two domains, (1,0,0) and (0,1,0) will be 
favoured. In the latter case, four nodes would be ob- 
served in the directions [±100] and [0 ± 10]. However, 
because the crystal symmetry of PrOs4Sbi2 is tctrahc- 
dral to begin with, any axial perturbation will lift the 
degeneracy of all three domains, any of which could be 
favoured. Therefore, in the A-phase, if all domains are 
present then tetrahedral symmetry with six nodes will 
be observed. Otherwise, only one domain is present, the 
symmetry will be D^C^), with two nodes. It is not likely 
that two out of three domains would be present in the 
A-phase, but could be possible if they were very close in 
energy. 

The same arguments also hold for the more compli- 
cated B-phase. Six domains are possible, with twenty- 
four non-degenerate nodes. If there is a single domain, 
then the symmetry is D 2 (E), and four nodes will be 
present. 



VI. CONCLUSIONS 

In this article, we have attempted to give a phys- 
ical description and comparison of the sequences of 
phase transitions normal-^ ^(C^) x K, — > D%(E) and 
normal— >D 2 (E), which we identify with the phase transi- 
tions seen in experiments, normal— > A— >B or normal— >B, 
respectively. Although this description is derived entirely 
from basic considerations of symmetry, a complicated gap 
structure emerges with several unusual features. First, 
the positions of the nodes in the B-phase are not located 
on any symmetry axes. Although this is allowed by sym- 
metry to occur in crystals with other point groups, such 
a feature has never before been considered. Second, be- 
cause the B-phase is triplet and non-unitary, there are 
two non-degenerate gaps. The only known example of 
this is Sr2Ru04, but in that case the two gaps remain 
close in energy4£ In PrOs4Sbi2, for a direct normal— >B 
transition, the energy difference is expected to be large. 
Finally, the proposed A^B transition, which is charac- 
terised by the splitting into two of the degenerate nodes 
of the the A-phase, is highly unusual. 

In summary, we have proposed phase transition se- 
quences in accordance with experimental evidence avail- 
able to date and studied its basic properties. Supercon- 
ductivity is best-described by a three component order 
parameter in the triplet channel. The superconducting 
phase has D 2 (E) symmetry, is non-unitary, and has four 
cusp nodes at unusual points on the Fermi surface. The 
presence of nodes leads to a quadratic dependence on en- 
ergy in the density of states, and power law behaviour in 
the specific heat and nuclear spin relaxation rate. There 
is also a second, higher energy, nodeless gap which may 
be experimentally accessible. 
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APPENDIX A: PROOF OF THE EXISTENCE OF 
NODES IN THE D 2 {E) PHASE IN THE TRIPLET 
CHANNEL 

In Section II, we found the gap function using ba- 
sis functions given by ([7]), and order parameter com- 
ponents (0, * I772 1 , |?7i|)- The gap function takes the form 
©, which vanishes either in the plane k y = at the 
points defined by \/\rii\ 2 b 2 — \r]2\ 2 a 2 k x = ±\rj 2 \bk z when 
1 771 1 2 6 2 > 1 772! 2 a 2 , or in the plane k z — at the points 
\JVl2Va? - \r]i\ 2 b 2 k x = ±\r] 2 \bk y when |r?i| 2 6 2 < |77 2 | 2 a 2 . 
In Section II, only p-wave pairing (basis functions linear 
in fe) was considered. In order to rigorously demonstrate 
the existence of nodes all possible higher order pairings 
must be included in the basis functions. We now consider 
this most general case. 

The most general form for the basis functions of the 
representation T in Th is 

d\ = (/ (k x , ky , k z ) , g(k x , ky , k z ) , h(k x , ky , k z )) (Al) 
d 2 = (h(ky,k z ,k x ),f(ky,k z ,k x ),g(ky,k z ,k x ))(A2) 

= (h'JW) 

d% = (g(k z , k x , ky), h(k z , k x , ky), f{k z , k x , ky) (A3) 

= (g",h"j") 

where /(fe) is odd in fe, g(k x ,k y ,k z ) is odd in k z and 
even in k x and k v , and h(k x , k y , k z ) is odd in k y and 
even in k x and k z . Eventually, we will find solutions to 
A_ (fe) = where one of the k's is zero (in agreement with 
the particular case of lowest order in k basis functions 
©), so we set /(fe) = now. 
Using dU), (©and (|AHIA3|) one finds 

A 2 _ = \vi\ 2 (g" 2 + h" 2 ) + \v2\ 2 (g' 2 + h' 2 ) 

- 2 1 771 1 1 772 1 \Jh" 2 g' 2 + g" 2 g' 2 + h" 2 h' 2 . (A4) 
Case 1: k y = 0: g" vanishes and 



A 2 _ = (\ m \h" - \mWg' 2 + h' 2 ) 2 . (A5) 

Nodes are found where A_ = 0, or where the function 

Mk x ,k z ) = h 2 (k z , k x , 0)-^U(g 2 (0, k z ,k x )+h 2 (0, k z ,k x )) 

(A6) 

vanishes. 

Case 2: k z — 0: h! vanishes and 

a 2 _ = (\mW - \mWg" 2 + h" 2 ) 2 . (A7) 
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Nodes are found where A_ = 0, or where the function 

4>2{k x ,k y ) = j^g 2 {k v ,0, k x )-(g 2 (0, k x ,k y )+h 2 (0, k x ,k y )) 

(A8) 

vanishes. 
We have 

Mk x ,0) = h 2 (0,k x ,0)-p^g 2 (0,0,k x ) 

mr 

&(0,fc«) = -pL h 2 {0,k z ,0)<0 
M 

h(kx,0) = P\^g 2 (0,0,k x )-h 2 (0,k x ,0) 

m\ 2 
= -<f>i(kx,0) 
<h(0,ky) = -g 2 (0,0,k y ) < 



If <fri(k x ,0) > 0, then <f>i[k x ,k z ) changes sign, i.e., there 
is a node of A_ in the k y — plane somewhere between 
the positions (fc x ,0,0) and (0, 0, k z ). Symmetry requires 
that there be (at least) four nodes on the Fermi surface. 
If ^>i(fca;,0) < 0, then 4>2(k x ,k y ) changes sign, i.e., there 
are four nodes in the k z = plane. 

Thus we have proved that, in general, the triplet phase 
with order parameter components (0, « [772 1 ? has four 
nodes in either the plane k y — or fc 2 = at the posi- 
tions [±a, 0, ±/3] or [±a, ±/3, 0], where a and j3 depend on 
the particular form of the basis functions. These nodes 
are "approximate" , in the sense that they are a conse- 
quence of symmetry and follow from the most general 
basis functions for the T representation. These nodes 
are also "rigorous", since the state (0, z 1 772 1 , couples 
to no secondary superconducting order parameters.— 
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